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ABSTRACT 

During hard X-ray observations of the Norma spiral arm region by the Nuclear Spectroscopic Telescope Array 
( NuSTAR ) in 2013 February, a new transient source, NuSTAR J163433-4738.7, was detected at a significance 
level of 8(7 in the 3-10 keV bandpass. The source is consistent with having a constant NuSTAR count rate over a 
period of 40 ks and is also detected simultaneously by Swift at lower significance. The source is not significantly 
detected by NuSTAR, Swift, or Chandra in the days before or weeks after the discovery of the transient, indicating 
that the strong X-ray activity lasted between ~0.5 and 1.5 days. Near-infrared imaging observations were carried 
out before and after the X-ray activity, but we are not able to identify the counterpart. The combined NuSTAR 
and Swift energy spectrum is consistent with a power law with a photon index of T = 4.1 4l 1 5 0 (90% confidence 
errors), a blackbody with kT = 1.2 ± 0.3 keV, or a Bremsstrahlung model with kT — 3.0 42 ' \ keV. The reduced-/ 2 
values for the three models are not significantly different, ranging from 1.23 to 1.44 for 8 degrees of freedom. The 
spectrum is strongly absorbed with A H = (2.8 42,3 4 ) x 10 23 cm~ 2 , (9 4l 7 5 ) x 10 22 cm' 2 , and (1.7 41 q 7 9 ) x 10 23 cm -2 , 
for the power-law, blackbody, and Bremsstrahlung models, respectively. Although the high column density could 
be due to material local to the source, it is consistent with absorption from interstellar material along the line of 
sight at a distance of 1 1 kpc, which would indicate an X-ray luminosity > 10 34 erg s . Although we do not reach a 
definitive determination of the nature of NuSTAR J163433-4738.7, we suggest that it may be an unusually bright 
active binary or a magnetar. 

Key words: Galaxy: stellar content - stars: variables: general - surveys - X-rays: individual (NuSTAR 
J 1634334738.7) - X-rays: stars 
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1. INTRODUCTION 

Hard X-ray surveys of the Galaxy provide an opportunity 

to discover populations of extreme sources. The promise of 

such surveys has been partially realized with the International 

Gamma-Ray Astrophysics Laboratory ( INTEGRAL ; Winkler 

et al. 2003), which carried out a 20-100 keV survey of the en- 
tire Galactic plane and has discovered hundreds of new sources 

(Bird et al. 2010; Krivonos et al. 2012), including new types of 
high-mass X-ray binaries (HMXBs), pulsar wind nebulae, and 

magnetic cataclysmic variables (CVs). The Nuclear Spectro- 
scopic Telescope Array ( NuSTAR ; Harrison et al. 2013), which 

launched in 2012 June and covers the 3-79 keV bandpass, is 

the first focusing hard X-ray telescope in orbit. While it has a 

much smaller field of view than INTEGRAL does, it has a much 
lower background and greatly improved sensitivity. Thus, one 

of NuSTAR ’ s science goals is to extend our view deeper into the 
Galactic plane to look for hidden hard X-ray populations. 


During its first year of operation, NuSTAR has initiated 
surveys of ~1 deg 2 areas both in the Galactic center (e.g., Mori 
et al. 2013; Nynka et al. 2013) and in a region that samples 
the spiral arm population in order to probe potentially different 
environments where X-ray binaries are found. A spiral arm 
region centered on Galactic coordinates of l — 33715 and 
b = 0° was chosen for having the highest known density of 
OB star associations (Russeil 2003) and HMXBs (Bodaghee 
et al. 2007, 2012). This is part of the Norma spiral arm region, 
which was identified early in the INTEGRAL mission as having 
an unusually high density of hard X-ray sources (Tomsick et al. 
2004; Dean et al. 2005; Lutovinov et al. 2005). The combination 
of active star formation and evidence that compact objects have 
already formed suggests that a survey by NuSTAR may uncover 
compact objects associated with populations of massive stars 
such as magnetars or faint HMXBs that are early in their 
evolutionary process and may have neutron star or black hole 
accretors. 
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Table 1 

X-Ray Observations and Count Rates for NuSTAR J1 63433-4738. 7 


Mission 

Instrument 

ObsID 

Start Time 

End Time 

Exposure 

Count Rate a 




(UT) 

(UT) 

(ks) 


Observations in 201 1 

Chandra 

ACIS-I 

12529 

Jun 16, 6.97 hr 

Jun 16, 12.53 hr 

19.0 

<0.19 

Chandra 

ACIS-I 

12532 

Jun 16, 23.85 hr 

Jun 17, 5.55 hr 

19.5 

<0.21 

Observations in 2013 

Swift 

XRT 

00080509001 

Feb 21, 21.70 hr 

Feb 21, 23.13 hr 

1.79 

i in+1.77 
A - 1U — 1.06 

NuSTAR 

FPMA 

40014004001 

Feb 22, 7.77 hr 

Feb 22, 17.52 hr 

19.4 

1.6 ±0.6 

” 

FPMB 

” 

” 

” 

” 

<1.3 

Swift 

XRT 

00080505001 

Feb 22, 11.88 to 

Feb 22, 13.71 hr 

1.92 

<0.54 

Swift 

XRT 

00080506001 

Feb 23, 0.83 hr 

Feb 23, 2.62 hr 

1.98 

<1.9 

NuSTAR 

FPMA 

40014007001 

Feb 23, 14.52 hr 

Feb 24, 1.77 hr 

22.6 

4.3 ± 0.6 

” 

FPMB 

” 

” 

” 

” 

4.4 ± 0.9 

Swift 

XRT 

00080508001 

Feb 23, 20.20 hr 

Feb 23, 21.70 hr 

1.97 

o o-fl.8 
Z.Z— i 2 

Swift 

XRT 

00080511001 

Feb 24, 13.72 hr 

Feb 24, 15.18 hr 

1.75 

<0.82 

Swift 

XRT 

00032728001 

Feb 28, 1.06 hr 

Feb 28, 6.24 hr 

4.95 

<0.45 

Swift 

XRT 

00032728002 

Mai- 3, 9.00 hr 

Mai- 3, 19.06 hr 

4.85 

<0.12 

Swift 

XRT 

00032728003 

Mar 5, 10.66 hr 

Mar 5, 20.75 hr 

4.47 

<1.1 

Chandra 

ACIS-S 

15625 

Mar 23, 8.30 hr 

Mar 23, 11.91 hr 

9.84 

<0.50 

NuSTAR 

FPMA 

30001012002 

Mar 23, 8.52 hr 

Mai- 23, 18.35 hr 

16.6 

<0.50 

” 

FPMB 

” 

” 

” 

” 

<2.0 


Note. a Count rates in the 3-10 keV band in counts per ks. The errors given are I rr and the upper limits are 90% confidence. 


The full surveys of the Galactic Center and the Norma region 
will be carried out over a period of ~2 yr; here, we report on 
the discovery of a transient source made during the first part of 
the survey. In this paper, Section 2 describes observations made 
with NuSTAR, Swift, and Chandra, as well as the procedures we 
used to reduce the data. The results are presented in Section 3, 
and we discuss possibilities for the nature of the transient in 
Section 4. 

2. OBSERVATIONS AND DATA REDUCTION 

As the first part of the NuSTAR survey of the Norma region, 
nine ~20 ks NuSTAR observations were performed between 
UT 2013 February 20 and 24. Each 13' x 13' field of view 
was partially overlapping with adjacent pointings, and the 
entire region covered was ~0.2 deg 2 . The results from all nine 
pointings will be reported in Bodaghee et al. (2014). Here, 
we focus on the observations that covered a new transient, 
NuSTAR J 163433 — 4738.7. These observations are listed in 
Table 1, including two that were obtained during the survey 
and a follow-up NuSTAR observation on 2013 March 23 that 
was coordinated with Chandra. 

In addition, ~2 ks Swift X-Ray Telescope observations of 
the region were carried out during 2013 February 21-24, and 
four of these observations covered NuSTAR J163433^1738.7. 
Table 1 lists these along with three other Swift observations 
that were acquired after our survey as part of another observing 
program. We also analyzed archival data covering the source, 
including two Chandra observations (ObsIDs 12529 and 12532 
with exposure times of 19.0 ks and 19.5 ks, respectively) that 
were acquired in 2011 as part of a survey of the same region 
being covered by NuSTAR (Fornasini et al. 2014). 

We reduced the NuSTAR and Swift data using HE AS OFT 
v6. 14 and the latest version of the Calibration Database files 
as of 2013 August 30. We produced cleaned event lists for the 
NuSTAR focal plane modules A and B (FPMA and FPMB) 
using nupipeline and for the Swift X-Ray Telescope using 


xrtpipeline, and further analysis of the event lists is subse- 
quently described. For Chandra, we processed the Advanced 
CCD Imaging Spectrometer (Garmire et al. 2003) data with the 
Chandra Interactive Analysis of Observations software, using 
chandra_repro to make event lists. 

We obtained near-infrared observations covering the 
NuSTAR J 163433-4738.7 error region. This includes J, H, 
and K s observations performed on 2011 July 19 with CTIO/ 
NEWFIRM in the framework of near-infrared mapping of the 
Chandra survey field. A detailed description of the data and 
their reduction can be found in F. Rahoui et al. (in prepara- 
tion). They were reduced with the dedicated IRAF package 
nfextern following the standard procedure — tailored for wide- 
held mosaics — which consists of bad pixel removal, dark sub- 
traction, linearity correction, hathelding, and median sky sub- 
traction. The resulting images were then flux-calibrated through 
relative photometry with the 2MASS catalogue. We also ob- 
tained A)- band imaging with the Ohio-State Infra-Red Imager/ 
Spectrometer (OSIRIS) at the Southern Astrophysical Research 
4.1 m Telescope. We used the f/7 camera, providing a 80" 
held of view, centered at the nominal coordinates of NuSTAR 
J163433— 4738.7. We obtained 36 x 60 s dithered exposures of 
the held on 20 1 3 April 3 under good conditions with seeing of 
0'.'7, and observed it for 9 x 60 s again on 2013 April 5, also 
under good conditions but with somewhat worse seeing of l'.'O. 
Reductions were done using the XDIMSUM IRAF package and 
photometric calibration was obtained by comparing to 2MASS 
All-Sky Point Source Catalog sources in the held. 

3. RESUFTS 

The new transient was discovered from an inspection of the 
image from the 22.6 ks NuSTAR observation that took place 
starting on 2013 February 23, 14.52 hr. As shown in Figure 1, 
the source was detected in FPMA and FPMB. To determine 
the signihcance of the detection, we extracted 3-10 keV counts 
from a 30" radius circle centered on the approximate position 
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Figure 1 . NuSTAR discovery images in the 3-20 keV energy band for NuSTAR J163433-4738.7. The source was detected in both of the NuSTAR focal plane modules 
A and B in a 22.6 ks exposure taken on 2013 February 23-24. The images have been rebinned so that the pixel size is 2 '! 5 and smoothed with a 6 pixel Gaussian. The 
scale on the bottom of the figure is in counts per pixel, and a logarithmic scaling is used. The apparent elongation of the source is due to the distorted point spread 
function shape at large off-axis angles. 

(A color version of this figure is available in the online journal.) 


of the source. We determined the background level using a 
nearby source-free circular region with a radius of 90". After 
background subtraction, we obtained 97 ± 14 and 99 ± 20 counts 
in FPMA and FPMB, respectively. The combined significance 
is 8. Ocr, confirming the detection. 

To determine the source position, we extracted all of the 
events from a 60" x 60" square region and made histograms 
of the counts, binning in the R.A. and deck directions. We 
performed x 2 fitting of the histograms with a model consisting 
of a constant (accounting for the flat background) and a Gaussian 
for the source. It should be noted that this is an approximation 
since the NuSTAR point spread function is non-Gaussian. We 
determined the centroids separately for FPMA and FPMB, and 
they are consistent with each other. The weighted average of 
the two centroids is R.A. = 16 h 34 m 33. s 42, deck = -47°38'4T.'9 
(J2000.0) with 3a statistical uncertainties of 6" 3 and 4"9 in R.A. 
and deck, respectively. After considering that the systematic 
pointing uncertainty for NuSTAR is ~8" (Harrison et al. 2013), 
the error region can be approximated with a 10" radius circle. 

We searched online catalogs (e.g., S1MBAD 16 ) for 
X-ray sources consistent with the position of NuSTAR 
J163433— 4738.7, but we did not find any likely candidates. 
The catalog from the 2011 Chandra survey (F. M. Fornasini 
et al., submitted) does not have any sources in the NuSTAR 
J163433-4738.7 error region. On the basis of this and the 
analysis subsequently described (including a reanalysis of 
the 2011 Chandra observations), we conclude that NuSTAR 
J 1 63433-4738.7 is a previously undetected source and that it is 
very likely to be a transient given the sensitivity of the Chandra 
observations. 

For the NuSTAR and Swift observations listed in Table 1, 
we estimated count rates or upper limits for NuSTAR 
J163433— 4738.7 using 30" radius source regions centered at 
the source position derived from ObsID 40014007001. The 
X-Ray Telescope 90% encircled energy radius is approximately 


16 http://simbad.u-strasbg.fr/simbad 


20", but we use 30" to account for the uncertainty in the source 
position as NuSTAR J163433-4738.7 is not bright enough in 
any of the Swift observations to improve the measurement of the 
source position. As described earlier, we used a larger, circular 
region that does not contain any detected point sources for back- 
ground. Selecting the background regions for NuSTAR requires 
some care because of scattered light from a nearby bright source 
(4U 163CM-7). 

The 3-10 keV count rates or limits obtained for all observa- 
tions are given in Table 1 . For NuSTAR, the background rates are 
high enough to use Gaussian statistics, and we use the Poisson 
limits tabulated in Gehrels (1986) for Swift and Chandra. In all 
cases, we quote the la error bars if the minimum of the la error 
region is positive. Otherwise, we give the 90% confidence upper 
limit. While the NuSTAR observation taken on February 23 pro- 
vides the only highly significant detection, the Swift observation 
with the highest count rate is the one that occurred during this 
NuSTAR observation. The only other possible evidence for ac- 
tivity from NuSTAR J163433-4738. 7 occurred on February 22, 
when NuSTAR obtained a 2.7a detection in focal plane module 
A; however, this is not confirmed by the focal plane module 
B data. 

For Chandra ObsIDs 12529 and 12532 (from 2011) and 
ObsID 15625 (from 2013), we analyzed the Advanced CCD 
Imaging Spectrometer data to search for a detection of NuSTAR 
J163433-4738.7. On the basis of an inspection, no sources are 
apparent in the 0.3-10 keV or 3-10 keV images. The source is 
7' and 5' from the Chandra aimpoint for the 2011 ObsIDs. At 
these off-axis angles, the 90% encircled energy fraction (EEF) 
radii (for 4.5 keV photons) are 7"4 and 4" 5 for ObsIDs 12529 
and 12532, respectively. ObsID 15625 was a dedicated pointing 
with the target on-axis, and the 90% EEF radius is 2". For 
ObsID 12529, the largest number of 3-10 keV counts within any 
7" 4 radius circle inside the NuSTAR error region is three, and, 
after accounting for background, we calculate a 90% confidence 
upper limit of <1.9 x 1 0 4 s -1 on the count rate (see Table 1). 
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Table 2 

Parameters for Fits to NitSTAR and Swift X-Ray Telescope Spectra 


14 


12 


» 10h 


10000 


Swift 


20000 
Time (s) 


30000 


40000 


Figure 2. 3-10 keV NuSTAR light curves (focal plane modules A and B 
after background subtraction) for NuSTAR J163433^4738.7 during ObsID 
40014007001. There is one data point per satellite orbit. The errors on the 
data points are at 1 a confidence. The time of a Swift observation is indicated. 
Zero on the time axis corresponds to MJD 56,346.60000. 


For ObsID 12532, the largest number of 3-10 keV counts within 
any 4" 5 radius circle inside the NuSTAR error region is two, and 
the count rate limit is <2.1 x 1 0 4 s _1 . For ObsID 15625, the 
largest number of 3-10 keV counts within a 2" radius circle is 
two, and the count rate limit is <5.0 x 1 0 4 s'. This is higher 
than for the 2011 ObsIDs because of the lower exposure time. 

We made NuSTAR light curves for ObsID 40014007001 
with several different time binnings between 0.1 s and 5500 s 
(the approximate satellite orbital period) in the 3-10 keV and 
3-79 keV bandpasses. The lack of any apparent variability in the 
0.1-10 s light curves rules out flares or bursts that might prove 
the presence of a neutron star. Figure 2 shows the 3-10 keV orbit- 
by-orbit light curve for FPMA and FPMB combined. At the 
5500 s binning, a/ 2 test shows that the 3-10 keV and 3-79 keV 
light curves are consistent with the source being constant 
over ~40 ks. 

Next, we extracted NuSTAR FPMA and FPMB spectra for 
ObsID 40014007001 and an X-Ray Telescope spectrum for 
ObsID 00080508001, which are the two observations with 
significant detections of NuSTAR J 1 63433-4738.7. These were 
fittedjointly by minimizing the Cash (or Q statistic (Cash 1979) 
using the XSPEC software package. The statistical quality of the 
spectrum is low, and it is well fit by a power-law, a blackbody, or 
a thermal Bremsstrahlung model. In all three cases, we included 
absorption using Wilms et al. (2000) abundances and Verner 
et al. (1996) cross-sections. The parameters are given in Table 2, 
and the errors quoted are 90% confidence for one parameter 
of interest, AC = 2.7. Although we used the C-statistic to 
determine the parameters, we also calculated the reduced-/ 2 
values, and they appear in Table 2. This quantity is slightly 
smaller for the power-law model (/ 2 = 1.23 for 8 degrees 
of freedom) compared with the blackbody model (/ 2 = 1.44 
for 8 dof) and the Bremsstrahlung model (/ 2 = 1.29 for 8 
dof). However, given the small number of dof, the difference 
in Xv is not significant, and the steeply falling power-law index 
(T = 4. C'iq) suggests that the emission probably has a thermal 
origin. The spectrum for the blackbody model is shown in 
Figure 3. 


Model 

A H a 

T or kT 

Flux’ 5 

C-statistic c 

Xv/ dof 

Power-law 

Blackbody 

Bremsstrahlung 

282l 3 4 d 

9 + _‘ 7 5 

172 1 , 7 

4 1 +I - 5 
—1.0 

1.2 ± 0.3 keV 

J o +2 - 1 
J -1 .2 

4+18 

^—2 

1 f) +0 ’ 8 
1,u -0.3 

1 6 +20 

9.5 

11.2 

9.8 

1.23/8 

1.44/8 

1.29/8 


Notes. 

a The column density in units of 10— cm 2 . 
b The 2-10 keV unabsorbed flux in units of 10 12 erg cm -2 s _1 
c We fitted the spectra by minimizing the Cash statistic. 
d All eiTors in this table are quoted at the 90% confidence level. 


We used the absorbed blackbody model and the count rates 
for NuSTAR , Swift , and Chandra given in Table 1 to calculate 
measurements or upper limits on the absorbed 3-10 keV 
flux, and these flux histories are shown in Figure 4. We also 
calculated the flux upper limits for the Chandra observations 
from 2011. The upper limits on the absorbed 3-10 keV fluxes 
are <1.6 x 10“ 14 ergcuT 2 s -1 and <2.5 x 10~ 14 ergcm~ 2 s _1 for 
ObsIDs 12529 and 12532, respectively. Although the count rate 
limits are considerably lower for these ObsIDs compared with 
Chandra ObsID 15625, the flux limits are similar because of the 
different effective areas for the ACIS-I and ACIS-S instruments. 
The lowest Chandra upper limit indicates that the flux from this 
source changed by at least a factor of 39^®, suggesting that it 
is a transient rather than simply being a highly variable X-ray 
source. 

The NuSTAR error circle includes dozens of near-infrared 
candidate counterparts brighter than K s ~ 1 7 (Vega magnitude 
system). In the NEWFIRM images from 2011, the brightest 
source in the error circle (2MASS J 16343288-4738393) has 
K s = 12.33 ± 0.05, H = 13.07 ± 0.06, and J = 14.61 ± 0.05, 
and the 2MASS magnitudes are consistent, indicating a possible 
lack of variability on long timescales. The next two brightest 
sources are 2MASS J16343362-4738479 at K s = 13.31 ± 
0.05 and the Spitzer/ GLIMPSE source G336.7870+00.01 1 1 at 
K s = 14.70 ±0.05. These may be more likely counterparts 
to NuSTAR J163433-4738. 7 because they are relatively highly 
reddened (H—K s = 2.35 ± 0.08 and 1.34 ± 0.08, respectively). 
However, we do not find any evidence that these sources are 
variable. For the two OSIRIS K s images taken on April 3 and 
April 5, we performed aperture photometry for all of the sources 
within the NuSTAR error circle, but did not find any variable 
sources. 

The high time-resolution of the NuSTAR data allows for a 
search for a coherent signal with periods P ^ 4 ms, covering 
the range expected for either an isolated rotation-powered pulsar, 
a binary, or a magnetar. Given the paucity of source counts in 
the observations listed in Table 1, we concentrate our attention 
on the ObsID 40014007001 data. Photon arrival times, adjusted 
for the NuSTAR clock drift, were corrected to the Solar System 
barycenter using JPL DE200 ephemeris and the NuSTAR derived 
source coordinates. We extracted photons in the 3-10 keV band 
from a 30" radius aperture centered on the source to optimize the 
signal-to-noise ratio. We searched for significant power from a 
coherent signal using an fast Fourier transform (FFT) sampled 
at the Nyquist frequency. The observation span was too short 
to consider the smearing of the pulse profile by spin-down of 
even the most highly energetic pulsar or the typical binary 
orbit period. The most significant signal found has a power 
of P = 35.18, corresponding to a probability of false detection 
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Figure 3. NuSTAR and Swift energy spectra for NuSTAR J 163433-^1738.7 on 2013 February 23-24 fitted with an absorbed blackbody model. The black, blue, and red 
data points and model lines correspond to the X-Ray Telescope, focal plane module A, and focal plane module B, respectively. The errors on the data points are 1 <j 
confidence. 

(A color version of this figure is available in the online journal.) 
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Figure 4. NuSTAR, Swift, and Chandra (marked with diamonds, filled cir- 
cles, and a square, respectively) absorbed 3-10 keV fluxes for NuSTAR 
J163433-4738.7 assuming an absorbed blackbody with a column density of 
Ah = 9 x 10 22 cm -2 and a temperature of 1.2 keV. The errors on the data 
points are la significance, and the upper limits are 90% confidence. 


of p = 0.8 for 2 25 FFT search trials. We conclude that no pulsed 
X-ray signal is detected in from NuSTAR J i 63433-4738.7. 
After taking into account the local background, we place an 
upper limit on the pulse fraction at the 3 a confidence level of 
f p < 36% for a blind search for a sinusoidal signal P > 4 ms. 

4. DISCUSSION 

In considering the nature of NuSTAR J 1 63433-4738.7, it is 
useful to estimate the X-ray luminosity that the source reached 
during its outburst. Although the distance is highly uncertain, 
the strong absorption indicates either a large distance through 


a region of the Galaxy with heavy extinction or absorption 
local to the source. The former is a strong possibility because 
the transient is in the direction of a region of the Galactic 
plane that is crowded with H u/molecular cloud regions. 
NuSTAR J163433-4738.7 is at Z = 336°787, b = +0?014, 
near a group of FI ii regions (336.732+0.072, 336.840+0.047, 
336.9-0.1, 337.147-0.181, and 337.3-0.1) at a distance of 
10.9 ± 0.2 kpc (Georgelin et al. 1996; Russeil 2003). 17 If 
the absorption is interstellar, then NuSTAR J 1 63433-4738.7 
is very likely beyond or in these molecular clouds, indicating a 
2-10 keV unabsorbed luminosity limit of >9 x 10 33 erg s' 1 
for the blackbody spectrum and >1.4 x 10 34 erg s 1 for 
Bremsstrahlung. In the following, we consider both possibilities: 
a distance of ~ 1 1 kpc and a peak luminosity of ~ 10 34 erg s _1 ; 
and a smaller distance with absorption local to the source and a 
lower luminosity. 

Considering the first possibility, a luminosity as high as 
1 0 34 erg s“ 1 would be unprecedented for at least two common 
types of X-ray sources in the Galaxy. Nonmagnetic CVs have 
strong optical outbursts along with persistent and transient X-ray 
emission. However, extensive studies have shown that their 
X-ray luminosities are in the 10 29-32 erg s~ 1 range (Baskill 
et al. 2005; Kuulkers et al. 2006). While magnetic CVs can 
reach higher luminosities (Kuulkers et al. 2006), they do not 
typically show outbursts. Second, active binaries, including RS 
CVn systems and low-mass flare stars, produce X-ray flares, 
some of which can last for about a day. Superflares that peak at 
X-ray luminosities of 1 0 32 ~ 33 erg s -1 have been seen 
(Franciosini et al. 2001; Osten et al. 2007, 2010; Pandey & 
Singh 2012), and there have been flares that have released 
>10 37 erg (Franciosini et al. 2001). IfNuSTAR J163433-4738.7 
is at ~1 1 kpc, then in addition to having a higher peak luminos- 
ity, the energy released is >4 x 10 38 erg, which is based on the 
source being at its peak luminosity for >40 ks (see Figure 2). 

17 This kinematic distance is estimated using the mean Galactic rotation 
curve, and the uncertainty does not account for possible deviations relative to 
the mean. 
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In summary, while we do not rule out the active binary possibil- 
ity (and see later in this paper further discussion on this topic), 
at the 1 1 kpc distance, the event seen by NuSTAR would need 
to be extreme for this explanation to be correct. 

A source type that might be a good match to the NuSTAR 
J 1 63433-4738.7 properties is the class of highly magnetic 
isolated neutron stars: magnetars. While these sources are best 
known for their very bright and brief (~0. 1 s) X-ray and gamma- 
ray flares, they also have persistent but variable emission that 
can easily reach 10 34 erg s' 1 or higher (Woods & Thompson 
2006). The X-ray spectrum of their persistent emission is 
dominated by a ~0.5-l keV blackbody (Woods & Thompson 
2006; Mori et al. 2013). Furthermore, as mentioned in Section 1, 
magnetars are associated with regions where high-mass star 
formation is occurring, such as the Norma region, and there is 
a known magnetar, SGR 1627-41, that is ~0’.'2 from NuSTAR 
J 1 63433^4-738.7. However, one possible counterargument to 
the magnetar hypothesis is that magnetar periods of activity 
usually last for months. From Figure 4, the peak activity from 
NuSTAR J163433-4738.7 could not have lasted for more than 
~1.5 days based on the Swift upper limits, and Chandra and 
NuSTAR place a tight upper limit on activity ~3 weeks after the 
NuSTAR detection. While we do not detect pulsations or ~0. 1 s 
flares, which would prove that the source is a magnetar, the 
pulsation search is limited by the statistical quality of the data, 
and flaring episodes are relatively rare. 

A blackbody spectrum could also be produced from an 
optically thick accretion disk in a black hole binary. Most 
often, such sources show temperatures of ~ 1 ke V in the inner 
parts of their accretion disks when they are at high luminosity 
> 10 36 — 37 er g s -i, which would require a very large distance 
of >100-300 kpc for NuSTAR 1163433-4738.7. Two sources 
(GRS 1758-258 and IE 1740.7-2942) have shown fainter soft 
states, but their blackbody spectra have been at temperatures 
of 0.4 keV (Smith et al. 2001) and 0.7 keV (del Santo et al. 
2005), which are lower than the 1.2 keV we see for NuSTAR 
J163433-4738.7. Also, transient black hole binaries typically 
have outbursts that last for several months, which is very 
different from NuSTAR J 163433-4738.7. 

Considering the second possibility that the source is closer 
than ~ 1 1 kpc, then the high column density must be due to 
material local to the source. INTEGRAL has found a large 
number of obscured HMXBs that have column densities of 
10 23- 24 cm 2 due to the wind from their supergiant companions 
(e.g., Walter et al. 2006). While many of the sources in this 
class are persistent, the supergiant fast X-ray transients (SFXTs), 
which have some members that are also obscured HMXBs, have 
outbursts that can be as short as a few hours, and it would not be 
unusual for an SFXT to produce X-ray emission that lasts for 
half a day (Smith et al. 1998; Negueruela et al. 2006; Romano 
et al. 2011). However, a possible problem for this interpretation 
is that the spectrum of NuSTAR J 163433-4738.7 is softer than 
usually seen for active SFXTs. During flares, IGR J17544-2619 
and IGR J 16479-45 14 show a power-law spectrum with a 
photon index near T = 1.4 (Rampy et al. 2009; Sidoli et al. 
2013). Somewhat softer spectra with T ~ 2. 3-2.4 can be seen 
during periods of weaker activity, but their spectra are only as 
soft as NuSTAR J163433-4738.7 when they are in quiescence. 

For CVs and active binaries, the challenge is to explain 
a column density as high as we measure with absorption 
local to the systems. One possibility is that material expelled 
from the system (either related to the evolutionary state of 
the stellar component or related to the X-ray flaring event) 


could cause increased column density. However, for one active 
binary (AX Ari), the column density was seen to increase to 
Ah = l.lxl O 20 cm' 2 during an event (Franciosini et al. 200 1 ), 
which is two to three orders of magnitude less than is required 
for NuSTAR J163433-4738.7. Variation in optical brightness 
of some stars has also been attributed to the ejection of material 
from the stars. An example of this is a change of 0.5 mag in 
the optical seen for KU Cyg (Tang et al. 2011). However, this 
would translate into a column density near 10 21 cm -2 , which is 
still far less than is required. 

The near-infrared information that we have provides only 
weak constraints on the source type. Given that we were not able 
to identify the counterpart, we can only say that the source must 
be fainter (before and after the X-ray flare) than the brightest 
source in the error circle (2MASS J16343288-4738393). This 
corresponds to K s > 12.3, H > 13.1, and J > 14.6, and these 
values are consistent with all of the possibilities for the source 
type discussed above. However, in the scenario where NuSTAR 
J163433-4738.7 is a nearby source with local absorption, these 
limits do provide some constraint. For example, a very nearby 
and luminous obscured HMXB is ruled out. 

In summary, NuSTAR J 163433-4738.7 is a fast X-ray tran- 
sient that has a thermal spectrum with relatively high tem- 
perature ikT = 1.2 keV for a blackbody or 3.0 keV for 
Bremsstrahlung). If its high column density is due to interstellar 
material, the source is probably distant (> 1 1 kpc), making the 
peak luminosity >10 34 erg s' 1 . We discuss the origin of the 
flare and suggest that the most likely possibilities if the source 
is distant are an unusually bright flare from an active binary or a 
short outburst from a magnetar. We also consider the possibility 
that the source is closer and that the absorption is local to the 
source. More NuSTAR observations in the Galactic plane will 
determine whether such transients are common and, hopefully, 
shed light on the nature of NuSTAR 1163433^1738.7. 

This work was supported under NASA Contract No. 
NNG08FD60C, and made use of data from the NuSTAR mis- 
sion, a project led by the California Institute of Technology, 
managed by the Jet Propulsion Laboratory, and funded by the 
National Aeronautics and Space Administration. The authors 
thank the NuSTAR Operations, Software, and Calibration teams 
for support with the execution and analysis of these observa- 
tions. This research has made use of the NuSTAR Data Anal- 
ysis Software (NuSTARDAS) jointly developed by the ASI 
Science Data Center (Italy) and the California Institute of Tech- 
nology (USA). R.J.A. was supported by Gemini-CONICYT 
grant 32120009. F.E.B. was supported by Basal-CATA 
PFB-06/2007 and CONICYT-Chile (through FONDECYT 
1101024, Gemini-CONICYT 32 1 20003 , and Anillo ACT 1101). 
L.N. wishes to acknowledge the Italian Space Agency (ASI) for 
financial support by ASI/INAF grant 1/037/12/0-011/13. The 
authors thank Harvey Tananbaum for providing Chandra Di- 
rector’s Discretionary Time for this project. This research has 
made use of the SIMBAD database, operated at CDS, Stras- 
bourg, France. 

REFERENCES 

Baskill, D. S„ Wheatley, P. J., & Osborne, J. P. 2005, MNRAS, 357, 626 
Bird, A. J., Bazzano, A., Bassani, L., et al. 2010, ApJS, 186, 1 
Bodaghee, A., Courvoisier, T. J.-L., Rodriguez, J., et al. 2007, A&A, 467, 585 
Bodaghee, A., Tomsick, J. A., Krivonos, R., et al. 2014, ApJ, submitted 
Bodaghee, A., Tomsick, J. A., Rodriguez, J., & James, J. B. 2012, ApJ, 

744, 108 


6 


The Astrophysical Journal, 785:4 (7pp), 2014 April 10 


Tomsick et al. 


Cash, W. 1979, ApJ, 228, 939 

Dean, A. J., Bazzano, A., Hill, A. B., et al. 2005, A&A, 443, 485 
del Santo, M., Bazzano, A., Zdziarski, A. A., et al. 2005, A&A, 433, 613 
Fornasini, F. M., Tomsick, J. A., Bodaghee, A., et al. 2014, ApJ, submitted 
Franciosini, E., Pallavicini, R., & Tagliaferri, G. 2001, A&A, 375, 196 
Garmire, G. R, Bautz, M. W., Ford, P. G., Nousek, J. A., & Ricker, G. R. 

2003, Proc. SPIE, 4851, 28 
Gehrels, N. 1986, ApJ, 303, 336 

Georgelin, Y. M., Russeil, D., Marcelin, M., et al. 1996, A&AS, 120, 41 
Harrison, F. A., Craig, W. W., Christensen, F. E., et al. 2013, ApJL, 770, 103 
Krivonos, R., Tsygankov, S., Lutovinov, A., et al. 2012, A&A, 545, A27 
Kuulkers, E., Norton, A., Schwope, A., & Warner, B. 2006, in Compact stellar 
X-ray Sources, ed. W. Lewin & M. van der Klis (Cambridge: Cambridge 
Univ. Press), 421 

Lutovinov, A., Revnivtsev, M., Gilfanov, M., et al. 2005, A&A, 444, 821 
Mori, K., Gotthelf, E. V., Zhang, S., et al. 2013, ApJL, 770, L23 
Negueruela, I., Smith, D. M., Reig, P., Chaty, S., & Torrejon, J. M. 2006, in The 
X-ray Universe 2005, ed. A. Wilson (ESA SP-604; Noordwijk: ESA), 165 
Nynka, M., Hailey, C. J., Mori, K., et al. 2013, ApJL, 778, L31 
Osten, R. A., Drake, S., Tueller, J., et al. 2007, ApJ, 654, 1052 


Osten, R. A., Godet, O., Drake, S., et al. 2010, ApJ, 721, 785 
Pandey, J. C., & Singh, K. P. 2012, MNRAS, 419, 1219 
Rampy, R. A., Smith, D. M., & Negueruela, I. 2009, ApJ, 707, 243 
Romano, P, La Parola, V., Vercellone, S., et al. 201 1, MNRAS, 410, 1825 
Russeil, D. 2003, A&A, 397, 133 

Sidoli, L., Esposito, P, Sguera, V., et al. 2013, MNRAS, 429, 2763 
Smith, D. M., Heindl, W. A., Markwardt, C. B., & Swank, J. H. 2001, ApJL, 
554, L41 

Smith, D. M., Main, D., Marshall, F., et al. 1998, ApJL, 501, L181 
Tang, S., Grindlay, J., Los, E., & Servillat, M. 2011, ApJ, 738, 7 
Tomsick, J. A., Lingenfelter, R., Corbel, S., Goldwurm, A., & Kaaret, P. 2004, in 
5th INTEGRAL Workshop on the INTEGRAL Universe (Noordwijk: ESA), 
413 

Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 
465, 487 

Walter, R., Zurita Heras, J., Bassani, L., et al. 2006, A&A, 453, 133 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 
Winkler, C., Courvoisier, T. J. L., Di Cocco, G., et al. 2003, A&A, 41 1, LI 
Woods, P. M., & Thompson, C. 2006, in Compact Stellar X-ray Sources, ed. W. 
Lewin & M. van der Klis (Cambridge: Cambridge Univ. Press), 547 


7 


